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ABSTRACT 

We determine the critical couplings for the deconfinement phase transition in 
SU(2) gauge theory on N T x A 3 lattices with N T = 8 and 16 and N a varying 
between 16 and 48. A comparison with string tension data shows scaling of the 
ratio T c j 1 ^fa in the entire coupling regime (3 = 2.30 — 2.75, while the individual 
quantities still exhibit large scaling violations. We find T c /yfa =0.69(2). We also 
discuss in detail the extrapolation of T c /Aj^g and -y/a/A^g to the continuum limit. 
Our result, which is consistent with the above ratio, is T c /Aj^g = 1.23(11) and 
V^/Aj^g = 1.79(12). We also comment upon corresponding results for SU(3) 
gauge theory and four flavour QCD. 
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1. Introduction 

Ever since the pioneering work of M. Creutz [1] the approach to asymptotic 
scaling, and thus the continuum limit, was one of the central issues in studies of 
gauge theories on the lattice. Although the first results were promising, it soon 
became clear that simulations on large lattices are needed in order to establish 
asymptotic scaling for asymptotically free quantum field theories such as QCD. In 
fact, recent numerical studies of the 0(3) cr-model in two dimensions [2] suggest 
that the asymptotic scaling regime may not be reached even for quite large cor- 
relation length while at the same time ratios of physical observables show scaling 
behaviour. 

The lack of asymptotic scaling as well as the scaling of certain ratios of physical 
observables has also been observed in SU(N) gauge theories. In particular in the 
case of the SU(3) gauge theory the deviations from asymptotic scaling are large 
and have been noticed early as a dip in the discrete /^-function [3,4], which led 
to deviations from asymptotic scaling by more than 50% for certain values of 
the gauge couplings. Although the dip is not that pronounced for SU(2) gauge 
theories, there are clear deviations from asymptotic scaling seen even for the largest 
values of the coupling, f3 = 2N/g 2 for colour group SU(N) : studied so far. In 
particular the analysis of the heavy quark potential and the string tension has 
been performed up to couplings as large as f3 = 2.85 [5]. Still there is no hint 
for asymptotic scaling at this large /3-value, which already corresponds to lattice 
spacings as small as ~ 0.05 fm. Furthermore, an analysis of the short distance 
part of the heavy quark potential suggests that the approach to the continuum 
limit may even be as slow as in the two-dimensional a-model [6]. 

One of the best studied quantities in SU(2) gauge theory is the finite temper- 
ature deconfinement phase transition. For lattices of size N T x N^. with A T < 6 
the critical coupling has been determined with high accuracy [7,8] and an extrap- 
olation to spatially infinite volume could be performed. Moreover, an analysis of 
the critical exponents at the transition point were in perfect agreement with those 
of the three-dimensional Ising model. Here it turned out that the Binder cumulant 
of the order parameter is an observable which is well suited to locate the critical 
coupling for given N T as finite spatial size corrections are only due to the presence 
of irrelevant operators. 

So far the analysis of the scaling of the ratio T c / ^fo was limited to a rather 
small coupling regime in the case of SU(2) as the critical couplings for the decon- 
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finement transition have been determined only for N T < 6. It is the purpose of 
this paper to further investigate the scaling properties of the SU(2) gauge theory. 
We will extend earlier studies of the deconfmement transition to lattices up to a 
temporal size iV r = 16. This will enable us to perform a quantitative test of scaling 
in SU(2). Furthermore, we can follow the apparent scaling violations over a large 
range of couplings, which allows us to analyze various extrapolation schemes to 
extract T c /A^g in the continuum limit. 

This paper is organized as follows: In the next section we discuss our strat- 
egy of calculating the critical couplings for the deconfmement transition on lattices 
with large temporal extent. In particular we discuss the finite-size scaling of cumu- 
lants of the Polyakov loop expectation value. In section 3 we present our numerical 
results. Section 4 is devoted to a discussion of scaling and asymptotic scaling and 
a detailed discussion of the extrapolation of these results to the continuum limit. 
Finally section 5 contains our conclusions. 

2. Finite-Size Scaling and the Continuum Limit 

Usually finite-size scaling (FSS) in the vicinity of a finite temperature phase 
transition is discussed for lattice SU (N) gauge models, without trying to make 
contact with the continuum limit, i.e., the scaling properties are studied on lat- 
tices of size N T x N d with fixed N T and varying N a , where d denotes the spatial 
dimension and the model is viewed as a d-dimensional spin system. In the con- 
tinuum limit the FSS properties of these non-abelian models should, of course, be 
discussed in terms of the physical volume, V = L d , and the temperature, T, in 
the vicinity of the deconfmement transition temperature T c . We will study here 
how the scaling behaviour of the continuum theory emerges from the lattice free 
energy on arbitrary lattices, i.e., when varying N T and N a . 

For a continuum theory having a simple critical point and a characteristic 
length L = V x l d the singular part of the free energy density, 

f=^ = - l ^l (1) 

is described by a universal finite-size scaling form [9,10], 

f s {T,H;L) = L- d Q fa (g T L^,g H L^ . (2) 

Here we assume that corrections to scaling from irrelevant scaling fields gtL Vi , 
proportional to negative powers of L, can be neglected; for N T = 6 a value of 



Hi = —0.9 has been found for the SU(2) gauge theory [8], showing that irrelevant 
contributions disappear rather fast with increasing N a . 

On a lattice of size N T x N% the length scale L and the temperature T are 
given in units of the lattice spacing 

L = N a a , (3) 
T -1 = N T a . (4) 

In general the lattice spacing a is a complicated function of the coupling (3 = 
2N/g 2 . The dependence of the lattice spacing on (3 is known only in the continuum 
limit in the form of the renormalization group equation 

Therefore it is advantageous to replace the length scale L by the dimensionless 
combination, 

TV 

L-T=— . (6) 

Using this ratio in the FSS relation for the singular part of the free energy density 
we get 



The scaling function Qf s depends on the temperature T and the external field 
strength h through thermal and magnetic scaling fields, 

g t = c t t(l + b t t) + 0(th,t 3 ) , (8) 
g h = c h h(l + b h t) + 0(th 2 ,t 2 h,h 2 ) , (9) 

with non-universal metric coefficients q, c^, bt and bh still carrying a possible iV T 
dependence. Here t is the reduced temperature, t = (T — T c )/T c , which in the 
neighbourhood of the transition point can be approximated by 



t=(p-p c ,oo) ' 



7_ Pc.oc 



(10) 



4iV6 

This approximation reproduces the correct reduced temperature in the contin- 
uum limit, which is easily verified by using Eq. (5). We note that t has 
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corrections to the leading term, which, however, contribute less than 8% in the 
relevant coupling regime, i.e., (3 > 2.0 for SU(2) and (3 > 5.5 for SU(3). This 
non-leading term introduces a logarithmic dependence of t, and thus Q/ s , on N T 
through Eqs. (4) and (5). A priori we cannot exclude that this violation of the 
otherwise universal N a /N T dependence is enhanced in the non-asymptotic scaling 
regime. However, as we will show later we do not find any hints for this in our 
Monte Carlo data. 

A non-vanishing magnetic field strength h corresponds to adding a symmetry 
breaking term of the form hZ(a,N T )N^P to the action. Here P denotes the 
Polyakov loop, defined as 



p=N- d j2 n u < 



(x ,x),0 



(11) 



x x = l 



and Z(a, N T ) is a renormalization factor necessary to remove divergent self-energy 
contributions to the Polyakov loop, which represents a static heavy quark source. 
A physical order parameter, < P p >, not vanishing in the continuum limit, a 
susceptibility Xp an d a normalized fourth cumulant (74 may then be defined through 
derivatives of f s with respect to the external magnetic field strength h at h = 0, 



< P p >= N T d Z(a,N T ) (P) = 



Xp 



df s 



dh 



7i=0 



94 = 




d 2 fs 




dh 2 


h=0 











(12) 



(13) 



(14) 



For the cumulant g\ the renormalization factors cancel and we end up with the 
usual expression known for the Binder cumulant of the order parameter [11,12,13]. 
The general form of the scaling relations derived from Eq. (7) is 



d n f s /dh-\ h=0 = ( ^ 



. 13+1 



(15) 



where the function Q n is defined as 



Qn(xi) = 



d r 



dx 2 ' 



:Qf 3 (xi,x 2 ) 



(16) 



x 2 =0 
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The finite-size scaling behaviour of higher order cumulants which are closely related 
to the derivatives defined in Eq. (15) is discussed in appendix A. We note that for 
<74, as well as the higher cumulants, the prefactors of the scaling functions cancel. 
Thus they take on unique fixed point values at gt = even on finite lattices, if we 
ignore corrections from irrelevant operators. The cumulants are thus well suited 
to determine the critical coupling from simulations on finite lattices. 

Let us first consider the case when N T is kept fixed. Then N T can be absorbed 
in the non-universal constants in g t and and we end up with the usual form of 
the finite-size scaling ansatz for f s . The critical coupling can be determined from 
the fixed point of g 4 (/3,N a ) [8,11,13,14]. 

Next we consider y = N a /N T fixed, varying N a and therefore N T accordingly 
as is needed to reach the continuum limit. Rescaling N a and N T by a factor b 
leads to a phenomenological renormalization (3((3;b;y) by the following identity 
for a scaling function Q 

Qi9MN T )(^y) = Qi^0 1 bN T )^y) . (17) 

Here we have kept explicit the dependence of gt on N T that comes in through the 
non-universal metric coefficients q and bt in Eq. (8). 

The property that the normalized fourth cumulant g 4 is directly a scaling 
function can be used to measure the discrete /3-function, A/3, by using the above 
identity for Q = g 4 . By writing Eq. (17) in terms of the scaling function Q and not 
for the scaling field gt directly we do not have to determine the metric coefficients 
in Eq. (8). The discrete /^-function is given by the shift in the coupling /?, which 
is necessary to get the same value of g 4 for the two different lattice sizes 

g 4 ((3-A(3 y ((3);N a ;N T )=g 4 ((3;bN a ;bN T ) . (18) 

The function Af3 y defined in Eq. (18) may depend on y through contributions from 
irrelevant scaling fields gi((3, N T )N^ 1 with negative exponents yi. The discrete (3- 
function can then be obtained by an extrapolation to y = oo 

A(3((3) = lim Ap y (p) . (19) 

The knowledge of (3 Cj00 (N T ) and A(3 allows to calculate the critical coupling for 
the rescaled lattice size 

l3 c ,oo{bN T ) = p Ct00 (N T ) + A(3 . (20) 
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In the following we will use this approach to determine /3 C;00 for N T = 8 and 16. In 
particular we will check the y-independence for iV r < 8 and use this information 
to justify our calculations for N T = 16 with moderate values for y. 

3. Critical couplings and the Binder Cumulant for N T = 8 and N T = 16 

As discussed in the previous section the Binder cumulants are well suited to 
determine the critical couplings for the deconfinement transition. Previously they 
have been used to determine the critical couplings for the SU(2) gauge theory on 
lattices with temporal extent N T = 4 and 6 [7,8]. We follow here the same strategy 
to determine the critical coupling for N T = 8. We have performed simulations on 
lattices of size 8 x N% with N a = 16, 24 and 32 at several values of (3 in the vicinity 
of the estimated critical point. We have used an overrelaxed heat-bath algorithm 
with 14 overrelaxation steps between subsequent "incomplete" heat-bath updates 
[15]. We used the Kennedy-Pendleton [16] algorithm with one trial per link. The 
acceptance rate was always larger than 96%. Most runs were done on the massively 
parallel CM-2 and so the use of a complete heat-bath algorithm would have led to 
a considerable waste of resources. Measurements were taken after each heat-bath 
update. Details on our parameter choices and number of iterations are given in 
Table 1, where we also give results for the estimated integrated autocorrelation 
times for the expectation value of the Polyakov loop. 

In Fig. 1 we show the Binder cumulant, which on a lattice is measured by 

(P 4 ) 

#4 = ^2 ~ 3 ' ( 21 ) 

with the Polyakov loop, P, given by Eq. (11). Also shown in this figure is an 
interpolation between results obtained at the various values of /?, which is based 
on the density of states method (DSM) [17]. Our implementation of the DSM takes 
into account each measurement and does not require the usual division of the range 
of the expectation value of the plaquette operator U p = N~ 1 TrU\U2Ulul in bins. 
This corresponds to an infinite number of bins and has the advantage to remove a 
free parameter, the number of bins, from the method. Because of the finite length 
of the Monte Carlo runs, the DSM provides reliable results only if histograms of 
Up belonging to adjacent couplings overlap. We convinced ourselves that at least 
2.5% of the data of one sample is contained in each of the overlapping tails of the 
distribution. 
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N a 


iV T 


P 


Nmeas 


Tint 


16 


8 


2.500 


25010 


3.4 


16 


8 


2.510 


21055 


3.7 


16 


8 


2.520 


21200 


3.6 


24 


8 


2.500 


55500 


10.6 


24 


8 


2.505 


29500 


10.0 


24 


8 


2.510 


29500 


13.5 


24 


8 


2.515 


39500 


9.8 


24 


8 


2.520 


14500 


14.9 


32 


8 


2.500 


30354 


29.6 


32 


8 


2.505 


34750 


17.1 


32 


8 


2.510 


38540 


15.3 


32 


8 


2.515 


28750 


18.0 


32 


8 


2.520 


18800 


30.0 


32 


16 


2.720 


30910 


10.9 


32 


16 


2.740 


40350 


12.6 


32 


16 


2.750 


35000 


12.0 


48 


16 


2.740 


43850 


45.0 



Table 1: Run parameters. 



The DSM allows an accurate determination of the intersection points of 
<74(/3; N a ; N T ) and g±(f3; bN a ; N T ), where N a is the spatial size of the smaller lattice 
and b is given by the ratio N a ' /N a . Taking also the largest irrelevant scaling field 
into account the coupling (3 c (N a ,b) at which two cumulants intersect varies with 
N a and b as [8,11,13] 

P c (N a ,b) = p Ci00 (l-ae) 

1/ (22) 

b 1 /" - 1 

Using the value v = 0.628 calculated for the three-dimensional Ising model [14] 
and yi = — 1 [8] we determine the critical coupling for iV T = 8 by an extrapolation 
to e = with the result 

Pc,oo( N r = 8) = 2.5115 ± 0.0040 . (23) 

The error has been obtained from a weighted linear fit to f3 c (N a ,b), where the 
error for the intersection points has been calculated using the jackknife method. 
The extrapolation is shown in Fig. 2. 
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Figure 1. The cumulant g 4 for N T = 8 and various values of the spatial 
lattice size as a function of the coupling (3. Solid curves are interpolation 
curves based on the density of states method. The dashed lines indicate 
the error on the curves estimated by the jackknife method. 

Taking into account the error band on the continuous curves obtained by the 
density of states method shown in Fig. 1 as well as the error on /3 C;00 we obtain 
for the fixed point value of g±, 

g 4 (N T ) = lim g^c^N^Nr) (24) 
g 4 (N T = 8) = -1.48 ±0.10 . (25) 

From the FSS theory as discussed in section 2 we expect, in fact, that g 4 is inde- 
pendent of N T , i.e., is universal also in the continuum limit. Indeed, this seems 
to be supported by our data. In Table 2 we give results for g 4 (N T ) obtained from 
simulations with N T =4, 6 and 8. 

All three values agree reasonably well within errors, although there seems to 
be a tendency to lower values of g 4 for larger iV T . At least partially this may 
be related to still too-low statistics for lattices with large A^ T as g 4 is a non-self- 
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Figure 2. The coupling value at the intersection points of g± for N a = 
16, 24 and 32 as a function of e which is defined in Eq. (22). The critical 
coupling fl c ,oo can be read from the figure as the section of the y axis at 
e = 0. 





9a 


4 


-1.38( 3) 


6 


-1.43( 8) 


8 


-1.48(10) 



Table 2: Fixed point value of g 4 for different iV T . 

averaging quantity [18]. We note, however, that a common value of ^ 4 taken from 
the spread of the data for N T = 4, 6, 8 and 16 shown in Fig. 3 at t(N a /N T ) 1 / 1/ = 
shows good agreement with the value of the three-dimensional Ising model [14,19] 

g A = -1.40(10) (26) 
Rising = -1.41(1) • (27) 

As discussed in section 2, g± is a scaling function, which in the continuum limit 
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depends only on the reduced temperature, t = (T — T c )/T c , and the dimensionless 
quantity y = LT. In terms of lattice variables, y is given by 

» = (28) 

and the reduced temperature in the asymptotic scaling regime is given by Eq. (10). 
For small values of t one expects, 

94(t, y) = gA{ty 1/v ) = g 4 + g±M lv . (29) 

This universal scaling behaviour indeed seems to be fulfilled quite well, as can be 
seen from Fig. 3, where we show g^{t,y) for various values of N T and N a . The 
critical exponent v has been taken from the three-dimensional Ising model [14]. 
Otherwise the presentation in Fig. 3 is parameter-free. We note that indeed the 
slope (74 i seems to be universal and within our accuracy does not show any further 
N T or even N a dependence. 

The universal scaling behaviour of g± in the vicinity of the critical point can 
be explored to determine the critical couplings for larger values of N T without 
an explicit determination of the fixed point g 4 from simulations on lattices of 
varying spatial size. We have calculated g& at three j3 values on a 16 x 32 3 lattice 
and at one f3 value on a 16 x 48 3 lattice. Details for these runs are also given 
in Table 1. In Fig. 4 we show the Binder cumulant as a function of f3. The 
critical coupling (3 C (N T = 16) can now be determined from the shift of (3 required 
to overlay the data of Fig. 4 with those for N a =16 and 24 shown in Fig. 1. A 
matching procedure according to Eq. (18) gives the shift values Af3 y=2 = 0.232(2) 
and Af3 y= 3 = 0.224(2). The two values show a noticeable difference, although 
the statistical error is relatively large. To investigate whether the observed y 
dependence is of significance, we checked if this is present also for N T =4,6 and 
8. We determined the shift Af3 y for y = 2, 3 and 4 using data from [7,8] and 
found that there is no significant y dependence of Af3 y . Furthermore Af3 y agrees 
with the shift A/3 = P c ,oo(N T ,i) — /3 c ,oo(N Tj 2) calculated from the infinite volume 
critical coupling. For the largest lattice size N T = 16 the critical slowing down 
is most severe, as is seen from T int given in Table 1. As mentioned before g± is 
a non-self-averaging quantity and, as for the susceptibility, the expectation value 
will be too small, if one has not enough independent measurements [18]. This 
seems to be true for our dataset on the N a = 48, A^ T = 16 lattice and we can not 
obtain the same precision as in the previous cases of A^ T = 4, 6 and 8. 
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Figure 3. The Binder cumulant #4 as a function ofty x l v for various lattice 
sizes as given in the figure. The critical exponent v has been taken to be 
the one of the three-dimensional Ising model, v = 0.628. For N T = 16 
and N a = 32 and 48 we also mark the error in x direction caused by the 
uncertainty in the critical coupling /3 C)00 . 

Averaging over the two values of A/3 y at y = 2 and 3 and assigning a large 
error, which includes both numbers, we get A[3(N T = 16) = 0.228(6), which is 
consistent with results obtained on lattices of size 32 [20]. Using Eqs. (20) and 
(23) this corresponds to a critical coupling, 

/3 c ,oo( N r = 16) = 2.7395 ± 0.0100 . (30) 

We want to point out that the combination of Monte Carlo simulation and 
the DSM is especially suited to determine the intersection points of #4 for small 
volumes. If the system becomes too large, the widths of the histograms of U p 
become very small and one needs a lot of simulations at different f3 values to cover 
a given range of the coupling. This was the case for our results on lattices with 
iV T = 16 and for these we could only use the single histogram version of the DSM. 
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Figure 4- The cumulant for N T and N a = 32 (dots) and 48 (star) as 
a function of (3. The straight line represents a linear fit to the data points 
for N a = 32. The four curves through the data points are obtained using 
the single histogram version of the DSM. The range of validity follows 
from the condition that at least 2.5% of the measurements in the data 
sample is contained in each tail of the distribution ofU p . This gives the 
range of validity in Up, which is converted to a coupling range using the 
computed energy coupling relation. The dashed lines indicate the error on 
these curves. 

4. Scaling and Asymptotic Scaling 

The critical couplings determined in the previous section can be used to test 
scaling and asymptotic scaling in SU(2) gauge theory over a much wider range 
of the coupling than was previously possible. Additionally we will compare the 
SU (2) results to existing data for SU (3) and four flavour QCD. After that we will 
try to point out common features of the data. 

Recently the heavy quark potential has been studied in the pure SU (2) gauge 
theory on quite large lattices and several values of the gauge coupling [5,6,21,22,23], 
which cover the range of couplings we have studied here. The string tension 
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extracted from the long distance part of the heavy quark potential can be used 
to test the scaling of dimensionless ratios of physical quantities. In particular we 
will discuss here the scaling of T c / yfa. In Table 3 and Table 4 we have collected 
results from most recent calculations of the string tension in SU(2) [5,6,21,22,23] 
as well as SU(S) [24,25] gauge theory. For four flavour QCD we used the value 
v^a = 0.332(2) [26]. 



N a 


iV r 


P 


\fo a 


8 


10 


2.20 


0.4690(100) 


10 


10 


2.30 


0.3690( 30) 


16 


16 


2.40 


0.2660( 20) 


32 


32 


2.50 


0.1905( 8) 


20 


20 


2.60 


0.1360( 40) 


32 


32 


2.70 


0.1015( 10) 


48 


56 


2.85 


0.0630( 30) 



Table 3: SU(2) string tension. 



N a 


iV T 


P 


y/a a 


16 


32 


5.60 


0.5295( 9) 


16 


24 


5.70 


0.4099(12) 


16 


24 


5.80 


0.3302(15) 


16 


24 


5.90 


0.2702(19) 


32 


32 


6.00 


0.2269(62) 


24 


32 


6.20 


0.1619(19) 


32 


32 


6.40 


0.1214(12) 


32 


32 


6.80 


0.0738(20) 



Table 4: SU(3) string tension. 



We have used a spline interpolation between these values to determine the 
string tension \faa((3) at the critical couplings of the deconfinement transition. 
These critical couplings are displayed in Table 5 and in Table 6 [27,28,29]. 
The resulting values for the ratio, 

^= = (V^a((3?(N T )) Mr)' 1 , (31) 
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Pc 


T C /A L 


^cM MS 




R 


2 


1.8800( 30) 


29.7(2) 


1.499(11) 


0.852( 6) 


3 


2.1768( 30) 


41.4(3) 


2.089(16) 


1.213(11) 


4 


2.2986( 6) 


42.1(1) 


2.125( 3) 


1.313( 3) 


5 


2.3726( 45) 


40.6(5) 


2.047(23) 


1.360(21) 


6 


2.4265( 30) 


38.7(3) 


1.954(15) 


1.354(13) 


8 


2.5115( 40) 


36.0(4) 


1.815(18) 


1.325(16) 


16 


2.7395(100) 


32.0(8) 


1.616(41) 


1.271(37) 



Table 5: Critical couplings for SU(2). 







Pc 


T C /A L 


T c/ A MS 




E 


12 


3 


5.5500(100) 


85.70 ±0.96 


2.975 ± 0.033 


1.217 ±0.038 


oo 


4 


5.6925( 2) 


75.41 ± 0.02 


2.618 ±0.001 


1.318 ±0.012 


oo 


6 


5.8941( 5) 


63.05 ± 0.04 


2.189 ±0.001 


1.299 ±0.002 


16 


8 


6.0010(250) 


53.34 ± 1.50 


1.851 ±0.052 


1.152 ±0.043 


16 


10 


6.1600( 70) 


51.05 ±0.40 


1.772 ±0.014 


1.155 ±0.012 


16 


12 


6.2680(120) 


48.05 ±0.65 


1.668 ±0.023 


1.110 ±0.017 


16 


14 


6.3830(100) 


46.90 ±0.53 


1.628 ±0.018 


1.111 ±0.016 


24 


16 


6.4500(500) 


46.27 ± 2.50 


1.537 ±0.087 


1.064 ±0.069 



Table 6: Critical couplings for 5*7(3). 



are shown in Fig. 5 as a function of aT c = 1/N T . Also shown in this figure is 
a result for four flavour QCD, which combines the measured critical coupling for 
N T = 8 (P c = 5.15 ± 0.05) [30] and the result for the string tension measured at 
P = 5.15 on a 24 x 16 3 lattice [26]. 

It is rather remarkable that scaling of T c /^/a is valid in the case of SU(2) to 
a high degree over the entire regime of lattice spacings from aT c < 0.25 (N T > 4) 
onwards. For 5*7(3) this seems to hold for aT c < 0.125 (N T > 8), although for 
larger values of the lattice spacing deviations are also only on the level of a few 
percent. Averaging over the region 0.0625 < aT c < 0.25 (4 < N T < 16) for the 
SU(2) data and over 0.0714 < aT c < 0.125 (8 < N T < 14) for 5*7(3) we obtain 



0.69 ±0.02 , SU(2) 

0.56 ±0.03 , SU(3) . (32) 

0.38 ±0.05 , QCD, n f = A 



A strong decrease of this ratio with increasing number of degrees of freedom (par 
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Figure 5. The critical temperature in units of the square root of the string 
tension versus the aT c = 1/N T for the case of SU(2) (circles) and SU(3) 
(triangles) pure gauge theory as well as QCD with four flavours of dynam- 
ical fermions of mass ma = 0.01 (cross). The straight lines correspond to 
one parameter fits to the data. 

tons) is obvious. Such a behaviour is expected, as a certain critical energy density 
can be reached already at a lower temperature when the number of partons is 
larger. Moreover, it is natural to assume that the string tension increases when 
the number of gluonic degrees of freedom increases. Both effects tend to lower the 
value of T c /y/a with increasing N and/or n/. 

A more natural scale to compare the critical temperature in pure SU (N) 
gauge theories is given by the glueball mass. The energy density in units of T 4 of 
a massive, free glueball gas is a function of x = mc/T only, 

3 oo „ 

e/T* = d^Y,[K 1 {nx) + -K 2 {nx)] , (33) 

n=l 

where d is a degeneracy factor, K\ and K 2 are modified Bessel functions. One 
thus might expect that the critical behaviour in the SU(2) and SU(3) theory is 
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controlled by the ratio T c /mc, with mg denoting the mass of the lightest glueball. 
In fact, using recent glueball data [21,31,32], we find that this ratio varies only 
slightly between SU(2) and SU(S), 

T c [ 0.180 ±0.016 , SU(2) . . 

m G \ 0.176 ±0.020 , SU(3) ' 1 ' 

It is well known that asymptotic scaling does not hold in the coupling regime 
considered by us, despite the fact that scaling works remarkably well. This suggests 
that there are universal scaling violating terms, which are common to the different 
physical observables and thus may partially be absorbed in a redefinition of the 
coupling constant. The need for such a resummation has been noticed quite early 
and various schemes have been suggested [33,34,35,36,37]. In the case of the SU(3) 
deconfinement transition a variant of the scheme originally introduced by Parisi 
[33] has been found to yield quite good results already for rather large lattice 
spacings [38]. Here the bare coupling, /?, is replaced by an effective coupling, (3e, 
which is related to the plaquette expectation value and thus takes into account 
rapid fluctuations in the pure gauge action at intermediate values of the coupling. 
Further details on the definition of (3e and a collection of plaquette expectation 
values needed to calculate (3e are given in appendix B. 

In Fig. 6 we show T c /A^g, using a conversion factor A^g/A^ [39] given by 

A, 



^ = 38.852704 exp {-^) • (35) 



In the case of full QCD with four flavours we used the value A^/A^ = 76.45 for 
the conversion factor. 

It is evident that the asymptotic scaling violations are much reduced when (3e 
is used as a coupling constant. This is particularly true for the case of SU(S) shown 
in Fig. 6b. The 0(g 7 ) term in the SU(N) /^-function would lead to violations of 
asymptotic scaling, which are of 0(l/lna). Our data for T c /Aj^g at finite values of 
the cut-off suggest, however, that the deviations from asymptotic scaling are well 
approximated by a correction term which is 0(a). This holds for both coupling 
schemes, with the bare coupling f3 and with the effective coupling /3 E . In the latter 
case, though, the coefficient of the correction term is much smaller. Performing a 
linear extrapolation of T c /A^jg to a = in the effective coupling scheme we find 

T c _ J 1.23 ±0.11 , SU(2) (36) 



A m \ 1.03 ±0.19 , SU(3) 
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Figure 6. The critical temperature in units of for the SU(2) (a) 
and SU(3) (b) gauge theory. Shown are data obtained by using the bare 
coupling constant (triangles) and the effective coupling (3e (circles) as in- 
put in the asymptotic renormalization group equation. The solid straight 
lines give linear extrapolations of these data sets to the continuum limit, 
a = 0, in the effective coupling scheme. The broken lines indicate a corre- 
sponding linear extrapolation using the bare coupling. The position of the 
filled symbols marks the result of the extrapolation. The dotted line in (a) 
marks how both coupling schemes approach in the continuum limit, where 
we assumed a linear form of T c / A-^ as a function of a given by the fit 
for the effective coupling scheme and used the third order expansion of the 
internal energy given in appendix B, Eq. (B.l). 
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where we used the region 0.0625 < aT c < 0.1667 (6 < N T < 16) for the SU(2) 
data and 0.0714 < aT c < 0.1250 (8 < N T < 14) for 5*7(3). We want to note 
that this result does not disagree with a linear extrapolation for the bare coupling 
scheme within errors. 

It is remarkable that T c in units of A^g differs only by less than 20% between 
SU(2) and SU(3) and also the corresponding value for QCD with four flavours, 
T c /A m = 1.05(20) [30], is surprisingly close to the SU(3) value. 

A corresponding linear extrapolation of the string tension y/a/Ay^ in the 
region 0.0625 < aT c < 0.1667 (6 < N T < 16) for SU(2) and 0.0714 < aT c < 0.1250 
(8 < N T < 14) for SU(3) to the continuum limit a = gives the result 



The data from Table 4 together with the extrapolated values are shown in Fig. 7. 
The corresponding value for four flavour QCD is yfa/Ay^ = 2.78 ± 0.50. This 
result is, of course, consistent with the scaling value of T c /y/a given in Eq. (32). 

We note that our value of the SU(2) string tension is only slightly larger than 
a recent estimate from the heavy quark potential [6], which yielded y/a/Ay^ = 
1.61(9). The agreement between these entirely different approaches gives addi- 
tional support for our a = extrapolation in the /^-scheme. 

5. Conclusions 

We determined the critical coupling for the SU(2) deconfinement transition 
on lattices with large temporal extent, N T = 8 and 16, using a FSS analysis of the 
Binder cumulant g^. We find 



A comparison with existing string tension data showed scaling in the entire cou- 
pling range 2.30 < f3 < 2.74 under consideration. For the SU{3) gauge group 
scaling sets in at a value of f3 ~ 6. We observe a 20% increase in T c / yfa when go- 
ing from SU(3) to SU{2), which can be understood in terms of the larger number 
of gluons in the SU{3) gauge theory, which tend to decrease T c and increase yfo 
at the same time. 




1.79 ±0.12 , SU(2) 
1.75 ±0.20 , SU(3) 



(37) 



/WiVr = 8) = 2.5115(40) , 
(3cM N r = 16) = 2.7395(100) . 
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In the bare coupling f3 = 2N/g 2 we observe no sign of asymptotic scaling 
up to (3 = 2.74. The alternative coupling scheme derived from the measured 
plaquette expectation value, shows much less deviation from asymptotic scaling. 
Our data indicate that these deviations are well described by O(a) corrections and 
a linear extrapolation of T c /Al and y/a/ A L to a = seems to be justified. 

We note, however, that our analysis clearly is not sensitive to 0(l/lna) correc- 
tion terms, which would be nearly constant in the rather small interval of couplings 
considered by us. 
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Appendix A: Cumulants of the order parameter 

We define cumulants of the Polyakov loop expectation value in the following 



way: 

(A.1) 



n t J dh n 



h=0 



On a finite lattice all odd cumulants are zero due to the Z(N) center symme- 
try. 

Ki = 0, i odd (A.2) 

K 2 =< P 2 > (A.3) 

K 4 =< P 4 > -3 < P 2 > 2 (A.4) 

K 6 =< P 6 > -15 < P 2 >< P 4 > +30 < P 2 > 3 (A.5) 
K 8 =< P 8 > -35 < P 4 > 2 -28 < P 2 >< P 6 > + 

420 < P 2 > 2 < P 4 > -630 < P 2 > 3 (A.6) 

Then moments of the order parameter can be expressed by cumulants: 

<P 2 > = K 2 (A.7) 

< P 4 > = K 4 + SK 2 (A.8) 

< P 6 > = K 6 + 15K 2 K 4 + 75K$ (A.9) 

< P 8 > = K 8 - 28K 6 K 2 - 35K% - 210K 4 K 2 - 1785K| (A.10) 

Now we consider three possible cumulant ratios given by: 
K 4 < P 4 > 

9 '=Kl = ^P^- 3 (A - 11} 

^ = f§ = ^W^s ~ ~ 75 (A.12) 

98 = it 4 = <P2 8 >4 + 28 ^ + 35^ - 210,74 - 315 (A.13) 
The general form of these cumulant ratios obtained from Eq. (A.l) is 

■ N x d(l-nm) 
77 

K n ~ ^.y Mi - n) (d n f s /dh-\ h=0 r ~ q 



K nm _ \777) d™f s /dh nm \ h=0 _ Q nm 



rn 
n 
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The ratio K nm /K™ is directly a scaling function. Therefore it can be ex- 
pressed as a function of since both are derived from f s and depend on the same 
argument g t (N a /N T ) 1 ' U . Then we expect g nm to be constant for a fixed value of 
9 a- 

In order to test the consistency of the FSS form for f s: Eq. (7), we measured 
the normalized moments M 6 =< P 6 > / < P 2 > 3 and M 8 =< P 8 > / < P 2 > 4 . 
From Eq. (A. 12) and (A. 13) we see that if g 6 and g 8 are constant then M 6 and 
M 8 should also be constant. The measured values of these moments together with 
the size of a typical error are given in Table 7. We see that M 6 and also M 8 are 
very stable for the different lattice sizes, thus supporting our FSS ansatz for the 
singular part of the free energy density. 



N a 


N T 


M 6 


M 8 


8 


4 


3.13 


6.94 


12 


4 


3.11 


6.82 


16 


8 


3.16 


7.17 


24 


8 


3.13 


6.97 


32 


16 


3.19(10) 


7.39(30) 



Table 7: Normalized moments at g± = —1.4. 



Appendix B: The effective coupling scheme 

In the Pe scheme the bare coupling, /?, is replaced by an effective coupling, 
Pe, which is related to the internal energy. The starting point is a perturbative 
weak coupling expansion of the internal energy, which is known up to 0((3~ 3 ) 
[40,41], * 

(u p ) = i - dp- 1 - C2 /r 2 - C3 /r 3 + o(/r 4 ) (b.i) 

Cl = (iV 2 -l)(l/4) 

c 2 = (N 2 - l)iV 2 (0.0204277- l/(32iV 2 )) 

c 3 = (AT 2 _ l)jV 4 (4/3) (0.0066599 - 0.020411/A^ 2 + 0.0343399/N 4 ) . 

* We thank H. Panagopoulos for informing us about the result of his 0(f3~ 3 ) 
calculation prior to publication. 
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An alternative coupling (3e can now be denned from the first two terms of (U p ), 



with Up denoting the plaquette operator as defined in the third section. The value 
of (Up) has to be determined by numerical simulations. Our data for the average 
action are listed in Table 8. For SU(3) the data of the average action were taken 
from large symmetric lattices [24,42,43,44]. They are summarized in Table 9. 



p 


i - (U P ) 


1.8800 


0.52637(5) 


2.1768 


0.43158(3) 


2.2986 


0.39746(1) 


2.3726 


0.37661(2) 


2.4265 


0.36352(1) 


2.5115 


0.34564(1) 


2.7395 


0.30869(2) 



Table 8: Expectation value of the average action for 577(2). 



p 


i - (Up) 


5.40 


0.52823(30) 


5.51 


0.50120(100) 


5.60 


0.47520( 2) 


5.70 


0.45100(80) 


5.75 


0.44105( 9) 


5.80 


0.43236( 5) 


5.90 


0.41825( 6) 


6.00 


0.40626( 2) 


6.10 


0.39592( 3) 


6.20 


0.38635( 1) 


6.30 


0.37788( 1) 


6.40 


0.36935( 1) 


6.60 


0.35438( 4) 


6.80 


0.34078( 1) 


8.00 


0.27935( 2) 



Table 9: Expectation value of the average action for SU(3). 
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The normalization in Eq. (B.2) is chosen such that the leading term in the 
weak coupling expansion of (3e = 2N/g E is identical to the bare coupling, 

Pe=P-c 2 /c 1 + 0([3- 1 ) . (B.3) 

Both couplings g and qe agree in the limit when g goes to zero 

9 2 e = 9 2 + (c 2 /c 1 )g 4 /(2N) + (cg/d)/ / \2Nf + 0(g 8 ) . (B.4) 

Additionally in both schemes the /9-function adg/da starts with the two universal 
coefficients 60 and bi, which can be seen from an expansion in the bare coupling g 

0^ = -b g 3 - b l9 5 + (3/2)(c 2 / Cl ) W + 0(g 7 ) 
da 

= -b g 3 E~b igE . (B.5) 

The constant C2/C1, however, redefines the lattice A-parameter in the /^-scheme 
relative to the /3-scheme, 

AT =eXp W • (R6) 
This yields a much larger A-parameter for the alternative coupling scheme, 

Ae _ \ 1.7217 , SU(2) 
A L ~ \ 2.0756 , SU(3) ' 

The /9-function in the /3e scheme can be separated into two terms 

a d9E =a dgd9E 
da da dg 



(B.7) 



dg E 1 16A^ 2 dE 



(B.9) 



dg g E g s N2-ld[3 

The first term is the usual /^-function while the second term dE/df3 is pro- 
portional to the "specific heat" c v . * The entire function dgs/dg asymptotically 
approaches unity as the coupling g goes to zero. 

For the SU(2) gauge theory a peak of the "specific heat" has been observed 
at (3 ~ 2.2 [45] which can be related to a nearby singularity in a generalized two 
parameter coupling space [46]. For SU(3) the same structure was found in the 
form of a huge bump in the region 5.2 < f3 < 5.8 [47]. 

As a consequence the dip in the discrete ,9-function of the bare coupling seems 
to be compensated by the peak of the "specific heat" , which motivates the choice 
of qe as a better behaved coupling. 



* Here we consider the system as a (d+ 1) dimensional model at "inverse tem- 
perature" P, which is not to be confused with the physical temperature T. 
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